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We evaluate the exact QED 2+1 effective energy for charged spin zero and spin 
half fields in the presence of a family of static magnetic field profiles localized 
in a strip of width A. The exact result yields an infinite set of relations 
between the terms in the derivative expansion of the effective energy for a 



^ ■ general magnetic field. Upon addition of the standard Maxwell magneto- 



static energy, the minimum energy configuration at fixed flux corresponds to 

a uniform magnetic field. 
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The effective action and the effective energy are important tools for the study of quantum 
electrodynamics [0,0. Using the proper-time technique, Schwinger showed [0] that the QED 
effective action can be computed exactly for either a constant or a plane wave electro- 
magnetic field. This result was later adapted to QED 2+1 with constant fields by Redlich || . 
For more general electro-magnetic fields one generally performs some sort of perturbative 
expansion, such as the derivative expansion Q which yields the large distance behavior of 
the theory. This approach was recently applied to QED 2+1 ||, showing (for example) for the 
special case of zero electric field and static magnetic field, that while the zero derivative term 
increases the effective energy, the next order correction term with two derivatives tends to 
decrease the effective energy. The question of vacuum stability is inaccessible in a derivative 
expansion, so more powerful tools are required to study it. In this Paper we make a first 
step towards the non-perturbative understanding of such a system by considering a new 
exactly solvable model which has a spatially-varying magnetic field. We expect our model 
to be relevant for recent investigations of symmetry breaking |],0] and finite temperature 
effects §0 in QED 2+1 . 

We show that the QED 2+1 effective energy for charged spin zero and spin half particles of 
arbitrary mass can be computed exactly in the presence of time-independent but spatially- 
varying magnetic fields of the form: 

B{x,y) = - — (1) 
[cosh(x/A)J^ 

Here, the magnetic field is localized in a strip of infinite extent in the y-direction and of width 
A in the x-direction. In the limit A — > oo the magnetic field in (0) tends to a uniform one 
of strength B. The constant strength B sets a length scale \ j\fe~B known as the magnetic 
length, and the derivative expansion regime corresponds to A >> 1/v ' eB. In practice, the 
system will be considered in a box in the y-direction of size L. The total flux $ of the 
magnetic field is then finite : $ = eBXL/n. We show that the effective energy has a simple 
exact integral representation involving elementary functions for all values of mass m, width 
A, electric charge e, and strength B. 

The motivation for the choice of these particular profiles for the magnetic field resides in 
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the fact that they may be fairly representative of the type of inhomogeneity that we expect 
in the system, while still being exactly solvable. They have been chosen here in the form 
of string-like flux tubes with finite width and infinite length. While we cannot, of course, 
analyze exactly the dynamics for arbitrary magnetic field configurations, the exact solution 
will permit us to study the system non-perturbatively within this family of profiles and gain 
rigorous constraints on the derivative expansion for more general fields. 

The starting point for the evaluation of the effective energy in the presence of the mag- 
netic field (1) is the functional determinant expression of the effective action, which we quote 
in Minkowski space-time : 

i J d?xC± = T In Det {D^ + m 2 + eEgF^ - it}. (2) 

Here, C± are the effective Lagrangians for bosons (+) and fermions (— ) of electric charge 
e in the presence of a gauge potential with = <9 M + ieA^. The Lagrangian £_ with 
Yf" = (i/4)[7 M , produces the effective action for a 4-component spinor consisting of 2- 
component spinors of masses m and — m respectively, whereas C + with = produces 
the effective action for spin complex scalars with mass m. Both systems are invariant 
under parity and time reversal. The determinants are understood to be regulated by Pauli- 
Villars masses in the ultra-violet, which will not be exhibited explicitly. Also, we calculate 
the effective action relative to that for zero electro-magnetic fields and therefore drop all 
contributions independent of the fields. The effective energy for static magnetic fields is 
then given by E± = — J dx dy C±. 

For the family of magnetic fields in (1), we may use the translation invariance of the 
problem in time and in the ^/-direction to work in an eigenbasis of frequency u and y- 
momentum k. Then the operator D^V 1 + m 2 + eYi± coincides with the Schrodinger 
operator of a solvable 1-dimensional quantum mechanical system. (Recall that the solvability 
of the constant field case is based on its relation to the solvable 1-dimensional harmonic 
oscillator system fTOfl ). To see this, we choose the gauge potential A x — and A y = 
\B tanh(x/A), which reproduces the magnetic field in (1). 

+ m 2 + eEgF^ = + V k (x) - u 2 (3a) 
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V k {x) = -js^l (tanh If] + \a 2 k {l + tanh f ) + \o?_ k {\ - tanh f ) (3b) 

Here, the following assignments for the parameters of this potential have been made, with 
a 3 denoting the spin projection eigenvalue +1 for spin up fermions and —1 for spin down 
fermions: 



a k = J(k - eBX) 2 + m? (4a) 



7± 



Wl + 4(eBA 2 ) 2 (+) bosons 

" (4b) 

| + eBX 2 a 3 (— ) fermions 
The 1-dimensional Schrodinger operator ([!a|) is of the supersymmetric quantum mechanical 
form [PH| , but with a fc-dependent super-potential Wk(x) = A-Btanh(x/A) — k. The /in- 
dependent potential Vk(x) is of the modified Poschl- Teller form, with asymptotic energy 
barrier heights — ► ±oo; the corresponding incoming and outgoing momenta are 
then just a±k, while u 2 is merely an overall shift in energy. The origin of the difference 
between the two asymptotic energy barriers may be simply understood in terms of classical 
electrodynamics. As a particle of charge e enters the magnetic strip from — oo, its kinetic 
energy is conserved while its momentum in the y direction behaves as p y (x) = p y (—oo) + 
e5A(tanh | + 1). If the momentum in the x direction is large enough, the particle traverses 
the strip, but for small x-momentum, it will instead be reflected off the strip and return to 
— oo. 

The evaluation of the determinants in (2) thus reduces to a spectral problem for the 
Schrodinger operator in (|3a|). Note that this Schrodinger operator has both a discrete and 
a continuous spectrum, in contrast to the constant B field case for which the spectrum 
is purely discrete. It is convenient at this point to analytically continue frequencies to 
imaginary values : uj — > iu, as usual; the effective energy is then given by 

r r+oo r+oo { d 2 \ 

£± = ±5? L dk L *" ln Det ("5? + ,4W + " ) (5) 

[Notice that we actually compute the difference between S± and the corresponding free field 
(B = 0) case; thus we may drop terms independent of B.\ Integrating by parts in uj and 
omitting ^-independent terms, the effective energy may be recast as 



2,L r+°° r+oo 

£± = =F— r / dk duo u 2 Tr G.^, k (6) 



AtT^ j — oo J — oo 

The resolvent Green function Gs,fc is the inverse of the Schrodinger operator, for general 
complex parameter E : 

(-^ + ^(*) - E ) J) = S(x - x') (7) 

It is standard to obtain GE,k(x,x') by matching the independent solutions of the homoge- 
neous equation, which are proportional to hyper-geometric functions in terms of the new 
variable £ = [1 + tanh(x/A)]/2 : 

Ul (x) = r (1 - ifF{z + \ - 7± , z + \ + 7 ±; 1 + 2a; £) (8a) 

u 2 (x) = r (1 - ifF{z + \ - 1± ,z + \ + 7± ; 1 + 2/3; 1 - (8b) 

Here, a and (5 are defined as the roots with positive real part of the following equations. 



a — f v a -k ~ E P = y®l-E z = a + [3 (9) 

With these conventions, u\ is regular at £ = (i.e. as x — > — oo), whereas «2 is regular at 
£ = 1 (i.e. as re — > +oo). The Wronskian of these solutions is a constant, given by 

w = *w.) - = l r( ^ ( j!gg:^ 7±) do) 

and the Green function G Ejk is therefore given by 

All information concerning the spectrum of the Schrodinger operator is contained in the 
trace of the Green function G Etk - We first compute 

/+oo j r+oo 

dx G Ejk (x,x) = — dx ui(x)u 2 (x) (12) 
-oo W J —oo 

Some care is needed in regularizing this integral at x — ±oo; for example, one can multiply 
the arguments £ and 1 — £ in the hyper-geometric functions in u\ and u 2 respectively by 
a factor 1 — e for e > and infinitesimal. With this regularization, the integral may be 
performed exactly in terms of the Euler psi-function ip(x) = T'(x)/T(x), and we find 



Tr Ge,u — — \ 



I I 

a /3, 



^{z + \- 7± ) + i>{z + \ + 7± )1 + f e (a) + f e {(3) (13) 



where / e (a) is a function of a, but not of f3, whose precise form is regulator dependent, 
but which does not contribute to the effective energy once we integrate over k. Here and 
in the following, it is understood that for the case of 4-spinors (— ) both spin states are 
to be summed over. The spectrum contains a finite number of bound states, which arise 
from the (simple) poles of the ^-functions in (]13) at z + | — 7 ± = — n for < n < 



7±- I/ 2 " y/m- a -k\ : 

E n = \{a 2 k + a\) - \-\n + \ - 7± ) 2 - ±\\a\ - a 2 _ k )\n + | - 7± )" 2 (14) 

The same discrete spectrum may be obtained by solving the homogeneous Schrodinger equa- 
tion (|7|) for real E and requiring normalizability of the eigenf unctions |12|]. The spectrum 
also contains a cut starting at o? k and another cut starting at at k , corresponding to the 
two barrier thresholds. For B = 0, the discrete spectrum is absent, whereas for constant 
magnetic field case (i.e. A — > oo), it reduces to E n = 2eB(n + 1/2) for bosons and 2eBn for 
fermions, as expected. 

We now complete the calculation of the effective energy, using the result of ( |T3| ) in the 
expression @ for the effective energy. First, since a only depends upon k + eBX, we may 
shift k by —eBX in the contribution of f e (at) in (|13|). Thus, the regulator-dependent f e terms 
in (p~3|) yield only 5-independent contributions to the effective energy in (||) and may be 
omitted. Next, we use the identity 

-Y{a + p) = dE (15) 
which suggests the change of variable from to to z = a + f3 and yields 



S+ = ±-4t- / <ik 



... dz-J (z 2 - (a k + a_ k ) 2 ) (z 2 - (a k - a- k ) 2 ) [^{z + \ - 7± ) + i[>(z + \ + 7± )] 

Z7T A J-oo J\a k +a^ k \ Z 

(16) 

The integration over k can be performed and we end up with the remarkably simple expres- 
sion 
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S ± = ±-^—l dz , Z{Z \ Z ° 2) = [^ + |- T± )+^+| + 7 ± )] (17) 



"4ttA 2 



z o ^z 2 - z 2 + (Am) 2 



where = XJ (eBX) 2 + m 2 . This expression can be rewritten in terms of elementary 



integrals by making use of the following representation of the ^-function [ 13| : 



ijj(x) = Inx - — - 2 [°° — 1 dt n r (18) 

^ V ; 2x Jo (t 2 + x 2 )(e 2nt -1) y 1 



The first two terms on the right hand side in ( |i~8| ) contribute -B-independent terms to S±. 

(It is necessary to sum over both spin states to see this in the fermion case.) For the third 

term, the z-integral can be carried out exactly and we obtain the following finite integral 

representation for the effective energy 

^ ,-r , L f°° , 1 . (A 2 m 2 + v V) , Am - iv± \ , , 

£ ± (L,mX,eBX 2 ) = —— dt — — \{b± - it y ^ln- — + c.c. 19 

47rA 2 Jo e 2nt ±l \ v± Am + iv ± J 

where c.c. denotes the complex conjugate, and 



(eBX 2 ) 2 + 1/4 (+) bosons 

(20a) 

eBX 2 (— ) fermions 



t 2 + 2itb + - 1/4 (+) bosons 

(20b) 

t 2 + 2itb_ (—) fermions 

Expression ( |T9"D gives the exact effective action for a background field (|XJ) and is the main 
result of this Paper. For definiteness, we now concentrate on the fermion case, but analogous 
discussions can be made for bosons. 

In the limit of vanishing mass, m = 0, the only relevant dimensionless parameter is eBX 2 , 
and so one finds an asymptotic expansion 

LX(eB) 3/2 ~ 1 r(j - 3/2)r(j + 5/2) 



in 



3=0 



(AneBX 2 y T(j + l)T(j/2 + l/4)r(3/4 - j/2) 



LX(eBf' 2 



C(3/2)-^C(5/2)J^ + ... 



(22) 
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where ((z) is the Riemann zeta function [13], The first term in this expansion agrees with 



the uniform B field case |3|Jl4]|, while the next term agrees with the first-order derivative 
expansion computation in IJ. 



For nonzero mass m, there is another dimensionless parameter, eB/m 2 , which is the 

ratio of the cyclotron energy to the rest mass energy. A double expansion of ([OJ) yields 

IJ B 2k+2j /2eB\ 2k+ > 

(2k)\ (2k + j-l/2)(2k + j-3/2)\m 2 ) 

(23) 



57T 



j=0- 



k=l 



with B n the n th Bernoulli number (JT3]] . Each power in A 2 corresponds to a fixed order in a 
derivative expansion of the effective action. The zeroth order term agrees with the eB/m 2 



expansion of the exact constant B field answer |j,|14], while the first order term agrees with 
the eB/m 2 expansion of the leading derivative expansion contribution found in ||. 

In fact, the specific configurations ([I]) give some insight in the more general case of a 
background magnetic field B(x) depending on one coordinate only. In a derivative expan- 
sion of the effective Lagrangian the terms with a total of 2j derivatives, have, up to 
integration by parts, a unique structure 



3 F [0]/eB^- 
\ ni 2 f 



£[2] =mi fl f £M£)' 



eB'(x) 



2j-2 



m 



3-2j p[2j] / eB(x) ■ 



eBW-Q\ 



x 



i=i 



m A 



eB'(x) 



m z 



1 i 



2j=4,6, 



(24a) 
(24b) 

(24c) 



where B®(x) denotes the Z-th derivative of B(x). Parity invariance forces Fi(x) to be 
even (resp. odd) for I odd (resp. even). A comparison with (|23D entirely determines the 
zero- and two- derivative terms (in agreement with ||) and gives at each order (2j > 2) a 
relation among the (2j — 2) functions f\ 2 ^\x). 

1 r(2fc-3/2) 



*n*) = -E 



^8vr 3 / 2 r(2fc + l) 
00 1 r(2A; — 1/2) 



k=l 



4tt 3 / 2 T(2k) 



B 2k (2x) 2k 
8 2k+2 (2x) 2k - 2 



1=1 s=l 



8-1 



„2Z+s-2 



r(3/2) 



dx 



s a 



3-1+ 



i + 2 (j-l + 2-s)a s r 



s-l 

l+l 



+ xa s j _ l+1 Fg j] (x) 



~ i r(2fc + j)r(2fc + j-3/2) 2k 

1^ o^ftr/ol, , T\ tVoTa "2fc+2j [2X) 



Z 8tt 3 / 2 T(2k + 1) 



T(2k) 



where a s p are the coefficients in the polynomial of degree 2p 
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, 1 2 P~ 1 V 

[l-t 2 )j- t = Y,a s p {-2) 2p - s ~ l t 2p - 2s (l-t 2 ) s (26) 

s =l 

As an example, the two functions appearing in the four- derivative terms (2j = 4) obey 

Thus, Eqs. (^) provide an infinite number of relations among the coefficients of the effective 
action for any background magnetic field with translation invariance in one direction. 

Finally, we may apply the results on the effective energy obtained here and in || to the 
study of the full QED2+1 theory In ||, the effective energy is obtained for large wavelength 
fluctuations in the fields. While the leading order term contributes positively to the energy, 
the next, two-derivative term contributes negatively. It was proposed in || that this behavior 
may drive the system towards a lowest energy state with inhomogeneous magnetic field. 
Indeed, since flux is conserved in QED 2 +i, it is natural to consider static fluctuations 5B(x) 
that leave the total flux unchanged. Under this constraint, the energy fluctuation is given 
by 

£{B + 8B{x)) - £{B) = X - J d 2 x {l + a (e 2 /V^B) [l - /3 (eBy^ddB/dB] 2 ] } (5B) 2 (28) 

where the coefficients a and (3 are positive functions of m 2 /(eB), except for bosons with 
sufficiently small mass as explained in ||. The coefficient (3 is typically of order 1, and 
within the derivative expansion regime (eB)~ 1 \d5B /5B\ 2 << 1, indicating that the constant 
background magnetic field is stable under variations that conserve flux. 

In the case of the special family of profiles for the magnetic field in (|l|) we can go beyond 
the leading order derivative expansion, using the exact effective energy ([T9|), together with 
the Maxwell term, to give the total energy as 

£* ot = LXB 2 + £±(L, mX, eBX 2 ) = + £±(L, mX, vr$A/L) (29) 

e 1 XL 

with flux $ = eBXL/ir. We have shown that the total energy ( f2"9"| ) at fixed flux $ is a 
positive monotonically decreasing function of the parameter A. As a result, a system with 
magnetic profile (|l|) is driven towards a system with uniform magnetic field. Thus, within 



this family, we have not found further support for our suggestion in || that magnetic field 
inhomogeneities lower the total energy at fixed non-zero flux. 

This raises a number of open questions, (a) Whether any other families of magnetic field 
profiles are integrable, within which non-uniform magnetic fields minimize the energy at fixed 
flux. The recursion relations Q2"5| ) suggest that magnetic field profiles with x-dependence only 
may be such candidates, (b) More generally, whether inhomogeneous magnetic fields can be 
found that minimize the energy for fixed flux. 

This work is supported in part by the NSF under contract PHY-92- 18990 and by the 
DOE under grant DE-FG02-92ER40716.00. 
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